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Abstract: We study the dynamics of fuzzy two-spheres in a matrix model which represents 
string theory in the presence of RR flux. We analyze the stability of the known static solutions 
of such a theory which contain commuting matrices and SU(2) representations. We find that 
the irreducible as well as the reducible representations are stable. Since the latter are of higher 
energy, this stability poses a puzzle. We resolve this puzzle by noting that the reducible repre- 
senations have marginal deformations corresponding to non-spherical deformations. We obtain 
new static solutions by turning on these marginal deformations. These solutions now have in- 
stability or tachyonic directions. We discuss condensation of these tachyons which corresponds 
to classical trajectories interpolating from multiple, small fuzzy spheres to a single, large sphere. 
We briefly discuss spatially independent configurations of a D3/D5 system described by the 
same matrix model which now possesses a supergravity dual. 
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1. Introduction 



Noncommutative spheres have made their appearance in string theory for a variety of reasons 
[4-12,15,16]. Apart from their intrinsic interest [15-22], non-commutative spaces could actually 
be a more realistic description of space-time at very short length scales [23,24]. One way to 
argue this would be that [12] in a curved space-time generically the NS B-field is present (by 
string equations of motion). The fuzzy geometries that we will encounter [4] arise because of 
the presence of RR fluxes through a three-sphere. One or the other such background field will 
generically be present in a curved background. Indeed, non-commutativity of space-time may 
be of a more basic nature, as indicated by its appearance in string field theories [32]. Fuzzy 
spheres have also been investigated in the context of dynamics of giant gravitons [8, 10,25-31]. 

In this paper we will study dynamics of fuzzy spheres, indeed more generic fuzzy two-branes. 
We will start with the matrix model action [33] with a Chern-Simon term, which arises due to 
background constant four form flux [4]. In section 2, we set up the model and its equations of 
motion. We then review how SU(2) representations with various spin, as well as commuting 
matrices, are solutions to the equation of motion. The nontrivial SU(2) representations geo- 
metrically correspond to fuzzy spheres. The energy of the various solutions suggests a picture 
of cascade, where reducible representations have higher energy than the irreducible represen- 
tation. In section 3, we study quadratic fluctuations around the solutions discussed in section 
2. Contrary to our expectations we find that all the solutions labelled by different SU(2) spin 
representations are all stable, i.e., no quadratic fluctuations have negative mass squared. This 
raises a puzzle regarding our cascade picture. We also briefly discuss some geometric aspects 
of the fluctuation spectrum. The puzzle about unstable solutions is resolved in an interesting 
way in section 4, by finding new solutions to the equations of motion, which corresponds to 
non-spherical configurations. ^ These solutions correspond to exactly marginal deformations 
of the original solutions. A subset of these solutions has earlier been discussed and analyzed 
in [35]. We study quadratic fluctuations around these new solutions and find tachyonic insta- 
bility. This is an indication that there exists at least one solution with lower energy. In the 
second part of Section 4 we discuss other new and interesting features of quadratic fluctuations 
around these new solutions. We show that in certain cases, due to deformation by a marginal 
parameter, number of flat directions increase. This jump in the dimensionality of moduli space 
is reminiscent of emergence of new marginal operators in c = 1 conformal field theory at self dual 
radius. We also discuss patterns of symmetry breaking in this section. In section 5, we relate 
our method with other approaches. We first compare our results with the spherical D2 branc of 
Bachas, Douglas and Schwcigert [6]. We also discuss Myers' dual D2 brane [4]. In section 6, we 
present the energy landscape more quantitatively and discuss dynamical evolution of multiple 
fuzzy spheres into a single big fuzzy sphere. In Section 7 we discuss effects of turning on a mass 
term in our matrix model. In section 7, we discuss relevance of our results to the dual SUGRA 
solution [9] and also to F1-NS5 system in type IIA string theory. 



^Morc details of deformations of fuzzy spheres and interesting issue of topology change due to non- 
commutativity will be discussed in detail elsewhere [34]. 



2 



2. The model 



The matrix model we will be concerned with is described by the action: 

S = To Tr[^Xi^ + ^[Xi,Xj][Xi,Xj] - '-k ei,kXi[Xj,Xk]) (2.1) 

where Xi,i = 1,2,3 are N x N matrices and Tq = \/1/k j is the zero-brane tension. Throughout 
this paper, we use units such that Itiol = 27r/^ = 1. This action must be supplemented by a 
Gauss law condition, arising from the Aq equation of motion. 

This model arises in several contexts. We will mention two cases (see Section 5 and Section 
8 for more details of these two and other cases): 

(a) Myers [4] discusses ( |2.lD to describe DO brane quantum mechanics (in type IIA theory) 
in the presence of a constant non-zero vev of the 4-form flux F^^jj^ = —Ineiju- Note that k 
has dimensions . The last term (the Chern-Simon term) in the action produces interaction 
between DO branes due to the 4-form flux. This interaction is demanded by consistency of the 
D-brane action with the T-duality symmetry of the string theory. 

(b) Alekseev et al. [12] derived this ( |2.1| ) as an effective action that reproduces the correlation 
functions of an SU(2) level k WZW BCFT, representing S'^ branes wrapping on an (with 
radius given by k). 

The approaches (a) and (b) are, in fact, connected. E.g., [6] shows that the WZW BCFT 
can be accurately described at large k in terms a DBI action for D2-branes in the presence a 
two-form RR flux and a gauge field background on the brane. These D2-branes are spherical and 
are embedded in ; as one takes the radius of the 3-sphere to infinity, the two background fluxes 
mentioned above exactly match with [4] description of those fluxes (in the D2-description). 

We will use this model as means to study the dynamics of non-commutative 2-branes, 
notably spherical branes. 

In the rest of this section we will review the static solutions of this action discussed in [4] 
and [12]. The static equations of motion are 

[Xj, {[Xi,Xj] - IK e.,fcXfc^] = 0. (2.2) 

This equation of motion admits an obvious solution 

[Xi,X,] = 0, (2.3) 

which represents N DO branes (with the coordinates given by the diagonal elements of the 
matrices Xi). In the absence of the Chern-Simon term it is well known that this configuration 
is a lowest energy configuration satisfying the BPS condition. However, as pointed out in 
[4], the Chern-Simon term modifies the situation radically. In particular, commuting matrices 
are no longer lowest energy configurations. In fact, among the available set of extrema, this 
configuration has the highest energy. 

Besides, it admits the following static solution which satisfies: 

[Xi,Xj]=iK eijkXk. (2.4) 
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Clearly any matrix representation of SU(2) will satisfy this equation of motion. It is easy to 
explicitly write down such Xj's 

Xi = Ui (2.5) 

where Ji,J2,J3 define, say, the dimensional irreducible representation of su(2). It is well- 
known that ( |2.5| ) define a fuzzy 5^, of radius r, where 

= = + 1), 2j + 1 = iV (2.6) 

Besides ( ^ ) and (|2.5D , [4], [12] also specify reducible solutions. That is, Xi can be a direct sum 
of several irreducible representations of SU(2). Such a configuration also solves the equation of 
motion. 

Xi = K ©^=1 4''^ (2.7) 



where 

s 

J2i2jr + l) = N (2.8) 

r=l 

It is clear from eq. ( |2.6| ) that this representation corresponds to s fuzzy spheres with radii 

R^r = ^^jr{jr + l) (2.9) 



In (2/7), the irreducible representations jr can also include the trivial representation Jj = 0. It 
is, therefore, obvious that ( |2.5D as well as (|2.3D are special cases of i \2.7\) . 

It is important to consider the energy of these static solutions which is given by the static 
hamiltonian, (in units of Tq, the DO-brane tension) 

V = Tr{-^[Xi,Xj][Xi,Xj] + eijkXi[Xj,Xk]). (2.10) 
The energy E of (|2.7[) is given by 



E = -ToK^ J ^ jrijr + l)(2jV. + 1). (2.11) 

The trivial representation jr = 0, Vr corresponds to the commuting set of matrices and the 
energy of this configuration is zero. It is worthwhile to mention here that we are measuring 
the energy of these configuration with respect to the energy of times the single DO brane 
mass (tension) . Clearly, zero energy for the commuting matrices is the reflection of the no force 
condition between BPS configurations. Nontrivial SU(2) configurations, however, have negative 
energy. This means these configurations are more stable than the trivial configuration. The 
lowest energy configuration is the spin j = (N — l)/2 dimensional irreducible representation. 

It is easy to illustrate this fact by taking simple low dimensional matrix examples. Let 
us consider an example of 4 x 4 matrices. In this case there are five distinct solutions to the 
equations of motion. The static solutions, their sizes, i.e., the radii of the fuzzy spheres and 
energies are summarized in the table 1 below. 



4 



Table 1 



Solution 


Radii of fuzzy spheres 


Energy 




(m units oi K/ivrj 


(m units OI Iqh ) 


(1) Commuting 


(0,0,0,0) 





(2) spin i © © 


(\/3/2,0,0) 


1 

4 


(3) spin i © i 


( V3/2, V3/2) 


-1/2 


(4) spin 1 © 


(V2,0) 


-1 


(5) spin 1 


Vlb/2 


-5/2 



2.1 Cascade 

As anticipated from the energy formula, the commuting matrices have highest energy and the 
irreducible representation has lowest energy. The intermediate solutions (2,3,4 in the above 
table) have energies in between these two extremes. These energy values suggest the picture of 
a "cascade", i.e., all the configurations except the irreducible representation should presumably 
have instabilities or "tachyonic" directions which should lead to their decay into the most stable 
configuration, viz. the irreducible representation. Geometrically this means the tachyonic insta- 
bility would set off the process of smaller spheres "fusing" into larger spheres which ultimately 
would turn into the largest sphere to achieve minimum energy. 

3. Analysis of instabilities 

Let us now look at the instabilities of these solutions. The quadratic fluctuation of H around a 
general static solution has the form 

5^H = i|- + To T¥( - [5Xi,X,][5Xi,X,] + [5Xi,Xi][5X,,Xj]) - 

2To Tr(([Xi, X,] - iKeijkXk)[5Xi, 5X,]) (3.1) 

= \['^'ia + yiaMij^abVib] (3.2) 

In the above we have parametrized the fluctuations 6Xi = J^a Via^a^i = 1, 2, 3; a = 1, 2, . . . 
— 1. Ilj are momenta conjugate to 5Xi while -Kia are momenta conjugate to yia- 
In order to understand the energy landscape of our model (cf. Sec |2.1| ) , it is important first 

to understand the neighbourhood of the critical points of the energy function. We will therefore 

be interested in The eigenvalues K^ujf^ of the quadratic fluctuation matrix Mij^ab (these are 

calculated for various n in Appendix A). 

We will find below that for all the critical points (|2.7| ) and ( |2.5| ) described so far, lo'^ > 0, 

hence there are no instabilities. This leads to a puzzle. We will describe first the calculation of 

eigenvalues and then return to the puzzle. 
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3.1 Spectrum of fluctuations 

For the solutions ( |2.7| ) or (|2.5| ), the second hne of ( |3.lD vanishes by virtue of equation of motion 
(|2.4| ). Using, further, the gauge fixing condition, [Xj,(5Xi] = 0, we get 

=i|-^[X„<5X/ (3.3) 

= \['^l + y^aMi%,y4 (3.4) 

The tables in Appendix A list eigenvalues in which marginal deformations Y , discussed in the 
next section, are turned on; in order to find the eigenvalues without them we have to turn them 
off. Let us mention, for example, the eigenvalue table for 3x3 matrices. We reproduce here the 
fluctuations around the spin 1/2 © and spin 1 solutions here. Absence of the y-deformations 
means that we have put c = 0: 

Table 2: TV = 3 
Numbers in square brackets refer to physical zero modes 



Solution 


00^ 


Multiplicity 







10 [3] 




2 


6 




(3/4) 


8 







8[0] 


1 


2 


6 




6 


10 



3.1.1 Zero modes 

For the spin 1 solution, corresponding to the largest sphere for Mat(3), the 8 zero modes are all 
of the form J^^^Xj = [Aa,Xj] where Aa,a = 1, ..,8 are SU{2>) Gell-Mann matrices. It is easy to 
see from the Gauss law = [Xj,nj] that these are gauge rotations. The notation 8[0] means 8 
zero modes, of which are physical. For the spin ^ ©0 the SU(3) rotations produce 7 gauge zero 
modes, since [Agj^j] = 0. There are three physical zero modes which are infinitesimal versions 
of the three y-deformations. It is easy to see the general formula for the degeneracy of exact 
zero modes is that there are a total of /Uq = — 1 + 2na zero modes out of which 2>na are 
physical (this excludes eigenvalues which become non-zero after y-deformations). 

3.1.2 Non-zero modes 

Around a spin j irrep ( |2.5| ), like the spin 1 example above, the eigenvalues k^w^ of the quadratic 
fluctuation M in ( |3.4| ) are given by 

ujf = ^l{l + l),l = l,...,2j (3.5) 

each with multiplicity = 2{l + 1). Thus, = 2 appears 6 times and = Q appears 10 
times for the spin 1 solution above. To see this spectrum, note that the quadratic fluctuation 
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operator in (3^) is then simply the laplacian (Ji)^ = {AdXi)"^. The Ji acts here by adjoint 
action on matrix-type fluctuations, rather than on column vectors, hence the eigenspaces split 
into representations I £ j x j = 0,1, 2j (we exclude / = in the above table since we have 
restricted to traceless 5Xi). The degeneracies /i/ are twice 21 + 1 since there are two independent 
fluctuations SXi,i = 1,2. 

By a simple generalization of the above argument, the eigenvalues around a reducible rep- 
resentation (e.g. spin j © j') include (a) the eigenvalue set around an irrep j, (b) the eigenvalue 
set around an irrep and (c) eigenvalues ujf = l{l -|- = \ j — + j' , with multiplicities 

//; = 4(2/ -|- 1). In the spin ^ © example above, we have oof^^ — 3/4 appearing = 8 times. 

3.1.3 Geometrical interpretation of eigenvalues 

To put all this more simply, the structure of the fluctuation, say 5Xi, around a solution spin 
j © f is 

5Xi= fj'"^' '^"^'^ (3.6) 



Recall that the solution is 



Xi = K( '^'^''^ ) (3.7) 

\ Oj'xj Jij'xj' J 

The eigenvalues coming from the block a, namely ujf = 1(1 + 1),^ = l,...,2j represent various 
multipole deformations of the fuzzy sphere carrying representation j. Similarly, block b rep- 
resents multipole deformations of the other sphere j' . The blocks c, represent deformations 
which involve both spheres and will be shortly identified with tachyonic directions which deform 
the spin j © j' solution to other solutions, e.g. the irrep j + j' ■ 

One can easily check that the results mentioned above apply to all the tables in the Appendix 

A. 



3.2 The puzzle 



We find that > around all the solutions ( |2.7| ) or (2.5) that we have studied so far. This 



is a puzzle. For one thing, it appears to contradict the cascade picture mentioned in Section 



2.1| . Furthermore it is also mathematically absurd to have a potential none of whose critical 
points have unstable eigenvalues. To put it differently, if all the solutions mentioned in table 
1 are locally stable extrema with no unstable direction then it points to existence of barriers 
separating different extrema labelled by SU(2) representations. If so then there has to be some 
unstable extremum separating two such stable solutions. It is clear from above analysis that 
if such an extremum exists it is not a SU(2) representation. One of the motivations for this 
work was to resolve this puzzle. As we will see in the next section, this puzzle is resolved in an 
interesting way. The resolution lies in the fact that there are new critical points which we will 
find. 
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4. New critical points 



In this section, we will reconsider the equation of motion ( p.2| ) and look for the general static 
solution and its properties. Notice that the equations of motion used so far, ( p.3p and (p^), for 
obtaining the solutions as SU(2) representations are themselves a kind of ansatz. It is easy to see 
that solution to these equations solve the original equation of motion. The converse, however, 
is not true. The original equation of motion supports many more solutions which, in general, 
are not SU(2) representations. Let us consider the following situation 

Xi = Xi + CiaYa (4.1) 

[Ya,Y^] = Q = [Ya,Xi]. (4.2) 



This is clearly a solution of the equation of motion ( |2.2| ) , whereas it does not solve ( |2.3D or (|2^ 



In the above, 1^, a = 1, 2, • • • , are any set of matrices satisfying (4^), and Cia are any 3na 
real numbers. Since we are concerned here with only traceless XiS we will consider only traceless 
y's. These constitute finite deformations of the solution Xi discussed in the last section. These 



finite deformations effected by 1^, by virtue of (4.2), do not change energy of the solution. 
Energy of this solution is the same as that obtained using Xi as a solution. The coefficients Cia 
give a continuous family of degenerate solutions, which in a sense parametrise the moduli space 
of these solutions. Another way of saying it is that the y-deformations are exactly marginal or 
flat directions. 

As an example, suppose Xi are given by the solution (3) of Table 1. We get a family of 
solutions labelled by c. The new solution (family) ([4. ID is, therefore 





Xi = K[ ' ~ \+CiY,Y =[ ^ \ (4.3) 



Here Jj = l/2cji are the spin 1/2 repn. of SU(2). In this example, there is only one marginal 
deformation Y . The specific choice Y = \ ® could be replaced hy Y = 1 ® diai which still 
satisfies (|4^). 

Now that we have the new solutions, let us now go back to the analysis of quadratic fluc- 
tuations (|3.1| ). As we will see, each of these new solutions represents a collection of spherical 
D-brane solutions of various sizes whose centres are separated by the y-parameters above. The 
particular case where the y-deformations refer to locations of DO-branes separated from a single 
spherical brane has been discussed in great detail earlier in [35]. 

4.1 Spectrum of fluctuations after y-deformation 



We now consider the eigenvalues of the quadratic fluctuation operator in (3T) after we turn on 
the y-deformations. Although the deformations do not change energy, the quadratic fluctuations 
around the deformed solution are rather different from those around the new solution. In 



particular, the third term in (3.1) which vanished by virtue of being proportional to the equation 



of motion (2.4) does not vanish any more. Earlier results on fluctuations can be obtained by 
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turning off Ya deformations. Quadratic fluctuations around a general solution are calculated in 
Appendix A. 

4.1.1 Tachyons 

As we can see, negative eigenvalues of 5^H appear around these new solutions. In some cases, 
the instabilities appear after a finite amount of marginal deformation; in some other cases, even 
infinitesimal deformations lead to solutions with tachyons. Let us consider the 4x4 matrix 
problem. In this case as mentioned in the previous section we have five inequivalent solutions. 
To show how tachyonic instability is obtained in different situations we will concentrate on three 
cases of spin 1/2 © © 0, 1/2 © 1/2 and spin 1 © 0. 

In case of single spin 1/2 representation, the background, i.e., the solution to the equation 
of motion (p.2|), is given by 




where, i = 1,2,3 and we have used Gellmann's notation for SU(4) generators. That is, Xi,i = 
1,2,3 are spin 1/2 generators (Pauli matrices) in the upper left 2x2 block (and zero in the 
rest). As and A15 are traceless Cartan generators which are proportional to identity in the 
upper diagonal 2x2 block. This solution geometrically represents a S"^ corresponding to spin 
1/2 representation and two isolated DO branes located at c and d. In the absence of marginal 
deformations, i.e., deformations along the flat directions, multiplicity of massless fluctuations is 
23. By turning on infinitesimal marginal deformation, four massless fluctuations acquire mass. 
While two of them have positive mass square, the other two are tachyonic, signifying instability 
of this configuration towards formation of more stable ones. These more stable configurations 
are obtained by allowing these tachyonic modes to condense. Instead of allowing these tachyon 
modes to condense, if we deform further along the marginal direction then for the deformation 
parameter {V3-V2) < 2c < (Vs + V^), we see that there is another tachyonic instability. One 
of the novel features of this instability is that due to the finite marginal deformation, some of 
the irrelevant operators become relevant or equivalently tachyonic. Later we will discuss other 
novel features of these new unstable directions. 

Now let us look at the spin 1/2 © 1/2 solution. This solution is represented as 



where, Cj denote marginal deformations. This is the first instance when one encounters a solu- 
tion consisting of two nontrivial SU(2) representations. This solution also develops tachyonic 
instability due to infinitesimal marginal perturbation. In this case also two marginal directions 
become relevant or tachyonic. Like the first case here also after a finite marginal deformation 
there is another tachyonic instability, where the marginal parameter the marginal parameter 
takes values in the interval 1 — ^1/2 < c < 1 + \/l/2. The situation is quite different for the 
solution 1 © 0. 

The spin 1 © solution is represented as 




Xi — Xi + CjAs + (ij.Ais 



(4.4) 



= (Tj © / + Cj/ © CJ3 



(4.5) 



Xi — Ji + CjAi5 



(4.6) 
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where, Ji consist spin 1 generators which occupy 3x3 upper left diagonal block and zero in the 
last diagonal entry. In this case there is no tachyonic instability due to infinitesimal deformation. 
Even more interesting feature is that no marginal deformations become massive. This fact 
is valid for finite marginal deformations as well. This configuration is not the lowest energy 
configuration and hence it is imperative to find the tachyonic directions in the moduli space of 
this solutions. As mentioned a moment ago, none of the marginal directions become tachyonic 
instead we find that irrelevant directions become relevant. Since there is no tachyonic instability 
due to infinitesimal marginal deformation and the first time we encounter the instability when 
the marginal deformation parameter 2c > \/6 — \/3, the stability radius r for this solution in the 
moduli space is given by r = (\/6 — \/3)/2. 

4.2 Enhanced symmetry 

Here we will comment on yet another novel feature of the fluctuation spectrum. As mentioned 
above every solution has a multi-dimensional moduli space. However, we will see later that some 
of these moduli are gauge artifacts. Remaining moduli parameters or zero modes are genuine, at 
least at the point in the moduli space where the solution can be written in terms of direct sum 
of SU(2) representations. In the previous subsection we saw that some of these zero modes do 
get lifted due to infinitesimal marginal deformations. We also encountered a situation when an 
irrelevant deformation became relevant after a finite amount of marginal deformation. In other 
words a deformation with positive mass squared becomes one with negative mass squared. A 
new massless mode which separates these two regions at a point in the moduli space produces 
a new marginal direction. This sudden jump in the dimensionality of the moduli space is quite 
reminiscent of emergence of new marginal operators in c=l CFT [36]. Consider an example of 
spin 1 (gi solution to the 4x4 matrix problem. This solution develops instability only when 
the marginal deformation parameter 2c > -v/6 — \/3- In fact, when 2c = \/6 — \/3 this solution 
has two additional marginal directions. 

4.3 Symmetry breaking patterns 

As indicated above, although the y-deformations do not change energy, they are not trivial 
symmetry directions. The table of eigenvalues in Appendix A follow the symmetry breaking 
pattern. Recall, in the absence of RR four form background, trivial SU(2) representation was 
the lowest energy (BPS) configuration. When N DO branes are located at one point in space, 
they give rise to U{N) local gauge symmetry. If these DO branes are moved away from each 
other, this U{N) symmetry breaks down to U{1)'^ . 

In the presence of RR four form flux, this picture changes dramatically. Now the trivial 
SU(2) representation has highest energy among possible extrema. Though it is a kind of maxima, 
there are no relevant directions when all the DO branes are coincident. It has a large number 
of zero modes as argued in the previous subsection, though several of them are gauge artifacts. 
Thus this extremum resembles, in the language of super Yang-Mills theory, the Coulomb branch, 
dimension of which is given by number of nontrivial zero modes. However, some of these flat 
directions are lifted as infinitesimal deformation along any of the nontrivial marginal directions 
reveals instability of this solution to formation of non-trivial SU(2) representations. Now let us 
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consider one nontrivial SU(2) representation, say, spin 1/2 and rest all singlets. This solution 
has lower energy than the trivial solution and in this case U (N) gauge symmetry breaks down 
to U{N — 2) X U{1) [9]. Like the trivial solution this one also has no relevant operators exactly 
at the point where we have SU(2) symmetry. Again we find that infinitesimal deformation along 
one of the several marginal (flat) directions show that there are relevant directions leading to 
more stable configurations. Generically, if we have one nontrivial representation of spin j and 
rest all are trivial representation then the gauge group is U{N — 2j — 1) x U{1). If instead we 
have several nontrivial representations then the residual symmetry is 



The U{1) symmetry associated with each fuzzy sphere gets enhanced if the reducible solution 
contains multiple number of SU(2) representation. For example, in case of 4 x 4 matrices, we 
get a solution which contains 2 SU(2) spin 1/2 representations. Apart from the U(l) symmetry 
of these configurations, we also have additional symmetry which rotates these two spin 1/2 
configurations into each other without any energy cost. This enhances the gauge symmetry 
to SU(2). The biggest sphere, which corresponds to j = {N — l)/2 spin solution is the most 
stable configuration. This solution has no nontrivial flat directions and in this case the gauge 
symmetry is completely Higgsed. Since there are no relevant as well as marginal directions, the 
spectrum has a gap and therefore this solution represents the Higgs branch. 

We will now illustrate this by taking an example of 4 x 4 matrices. Among the five classical 
solutions, the trivial solution preserves C/(4) gauge symmetry. The solution with spin 1/2 ©0©0 
has U{2) <Si U{1) symmetry, where U{2) is the symmetry of two coincident DO branes and U{1) 
is contributed by the fuzzy sphere associated with spin 1/2 representation. The extremum 
corresponding to spin 1 © breaks the symmetry down to U{1) © U{1). The minimum energy 
solution corresponds to spin 3/2 and in this case there is no residual gauge symmetry left over 
in the problem. As mentioned earlier, 4x4 matrix also has a solution with spin 1/2 © 1/2. 
Since the spin 1 /2 representation appears twice in this solution we have an enhanced symmetry, 
which rotates two spin 1/2 representations into each other. This results in an additional U(2) 
gauge symmetry of the solution. 

5. Comparison with other Approaches 

In this section we will mention other related works and compare the results obtained here with 
those obtained in the other approaches. For facility of comparison, we will reproduce our results 
from Sections 2 and 3 about the radius Ri\f and the energy of the irrep and the eigenvalue 
wf of the quadratic fiuctuation operator M (|3.l|) : 



m 



U{N-Y,{2ji + l))^U{l) 



= -^^'NiN^ - 1) 

= i^'^'^i = + 1)> / = 0, 1, . . . , iV - 1 



(5.1) 



(5.2) 
(5.3) 
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5.1 Spherical D2-branes in S3 x 

This model is studied in [6]. Strictly speaking, this is not a valid background of string theory 
(has dilaton tadpoles etc, for instance), but for the purposes of the present section this fact can 
be ignored. The string theory is described by a level-A; SU(2) WZW CFT, (the D2-branes are 
described by boundary states of this CFT). [6] also presents an alternative discussion in terms 
of a Born-Infeld-Chern-Simons (BICS) action for D2 brane. We will discuss the latter first. 

5.1.1 Born-Infeld (BICS) description 

This description uses the following closed string background 

ds'^ = —{dip^ + sin^ 7p{de^ + sin^ edcp"^)) 

k , sin2w!', 

B = —{^ sin edOdcj) 

2tt 2 

N 

F = -—sin eded(l) (5.4) 

The region of validity of the DEI description is [6] : 

A;> l,iV> 1,A;/A^> 1 (5.5) 

As [6]shows, a spherical D2-brane in such a background, wrapping some S'^ G is energetically 
stable, provided that the area of the is chosen as 

k TT 

Area = 47r— sin^ = —r- (5.6) 
zvr k 

Here A^ refers to the quantized i^-flux, or equivalently, the number of DO-branes. This corre- 
sponds to a radius 

In the region ( |5.5| ), this agrees with ( |5.1| ), provided we choose 

Note that we are working in the convention 

27ra' = 1 (5.9) 
The mass of the D2-brane turns out to be [6] 

Mn = NTo + En, 

En = 2kT2sin^ - NTo = -To^^(l + 0{jf) (5.10) 



This also agrees with (|5.2D under (5^) and (|5.5|) . Note that in our convention (|5.9D Tq 
4TTa'T2 = 2ttT2. 

The eigenvalues of quadratic fluctuation operator are 



^?,^ = ^^^,^ = 0,l,...,Ar-l (5.11) 



which agrees with (5.3). 
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5.1.2 BCFT description 

The BCFT results for the energy and the quadratic fluctuation eigenvalues are [6] 

Nit 

En = 2{k + 2)T2sm{^-^) (5.12) 

-3±}L (5 13) 

It has already been shown in [6] that these results agree with the Born-Infeld results in the region 

5.1.3 The ARS matrix model 

The ARS matrix model [12] is related to the formulation of the above system as a bound state 
of DO-branes. The matrix model is given by 

S = ToTri^^j + l[cl>„ cl>^f - ^e^i^$,$^.<I>,)] (5.14) 



This matrix model itself agrees with the matrix model that we have presented ( |2.lD , provided 



relate $j to the Aj as follows 



/ k [2/K 
$. = ^-A. . = ^- (5.15) 

The results for the radii, energy and the fluctuation spectrum of course should agree since the 



matrix models themselves agree. The relation between k and k is the same as in (5.( 
5.2 D2-branes in II A in the presence of RR flux 

This is the situation considered in [4] . The matrix model (DO-description) is identical to the one 
in ( p.ll) . The dual formulation in terms of the world- volume theory of the D2-brane uses BICS 
action in the presence of a RR-flux background 

Ftijk = -2f€ijk (5-16) 

and a U(l) flux on the world- volume 

N 

F = —sinede Ad(j) (5.17) 



We should identify 

The BICS action is given by 



f = K (5.18) 



As discussed in Myers, this action has extrema when 



Nf Nk , , 
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At these values of the radius of the D2-spheres, the energy evaluates to be 



EN = -^N'f (5.21) 



Under the identification ( ^.18 ), we get perfect agreement with the results of our matrix model 
to leading order in 1/N. 

Regarding the fluctuation spectrum, it is easy to work out the quadratic fluctuations of the 
BICS action around the spherical solution for the D2-brane. The calculation is similar to the 
one presented in [6]. The result is 

^ - R 2R ^^-^^^ 

Thus, the fluctuation spectrum is of the form of a Kaluza-Klein spectrum on a sphere of effective 
radius R' . Because of the relation between R' and R, the KK spectrum goes as R in stead of 
1/R. The fact that the KK spectrum can be modified because of noncommutativity has been 
found elsewhere, e.g. in [28,30,37]. Because of ( |5.20| ), the above fluctuation spectrum ( |5.22| ) 
exactly agrees with our result (157 



6. Dynamics 

In the previous sections we developed a fair idea of the potential energy landscape of our model. 



We found in Section 3 that the cascade picture as suggested in Section is not quite valid; 
however there are exactly marginal deformations of the reducible representations (Section 4). 
In terms of the energy landscape these flat directions are like ridges; we found in Section 4 that 
if we start from any critical point corresponding to a reducible representation, and move along 
the ridges, then downward slopes develop along which one can roll down to lower energy critical 
points. Clearly such a landscape offers rich dynamics. In the present paper we will not attempt 
to solve the full dynamics of the X^{t) which is a formidable problem; we will instead focus on 
small submanifolds of the configuration space which captures important features of the above 
landscape and discuss the dynamics within such ansatzes. 

To be concrete, let us consider dynamical evolution from one SU(2)-representation Jj an- 
other, say Jj. We will discuss this dynamics within the following ansatzes in turn. 

6.1 One-parameter ansatz 

We wish to solve the ( p. 2]) with the boundary condition 

Xiito) = Ui (6.1) 
Let us try the following ansatz (c/. Bachas-Hoppe-Pioline [38]) 

Xi{t) = K{ji + f{t){Ji-J,)) (6.2) 
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where Jj, Ji are two SU(2) representations which we wish to connect by a dynamical trajectory 
parametrized by f{t), ( |6.lD translating to 



/(to)=0,/(ti) = l. 



(6.3) 



We note that for Jj = 0, the above ansatz is similar to the spherical ansatz of the giant graviton 
scenario, where the energy is viewed simply as a function of a single parameter of the (DO or 
D2) configuration, namely the radius. Since the flat directions and the resulting new solutions 
of Section 4 represent non-spherical deformations, we will need at least a two-parameter ansatz 
to reproduce the physics of the flat directions. 

To begin with, however, we will stick with the simple ansatz ( |6.2| ). We discuss the case where 
Ji is the N = 2j + 1-dimensional irrep, with even, and Jj corresponds to spin j' (B = 
{j -l/2)/2,i.e.N' = 2j' + l = N/2. 

Using (K^ in (O) leads to 



/ = -V'if), Vj{x) = -aj + hjx^ - Cjx^ + djx^ 



(6.4) 



where the coefficients depend on the spin j. aj, bj,Cj and dj are found to be positive. 

The above problem (S.4) is that of an asymmetric double well potential (see Figure (|l|)) if 
we consider f(t) as the position of a particle on a line. 




Figure 1: The asymmetric double-well for j = 3/2, 5/2, 7/2 

In Fig. (H), the representation Jj, corresponding to the double-sphere, is the local minimum 
on the left whereas the irreducible representation Jj corresponds to the absolute minimum. Note 
that within this ansatz it is not possible to classically evolve from Jj to J,. This is of course 
similar to the puzzle referred to in Section 3. Just as that puzzle was solved by the discovery of 
flat directions, we need here to include one more parameter in our ansatz. 

6.2 Two-parameter ansatz including marginal deformations 

The asymmetric double-well potential, arising from the ansatz ( |6.2D misses the flat direction 
present in the actual problem. Let us try a new ansatz which has a new term involving the flat 
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direction Y (cf. Q) 

Xi{t) = + f{t){Ji - Ji) + g{t)d^) (6.5) 
subject to the boundary condition 

f{to) = 9{to) = 

/(ti) = = (6.6) 

The action ( |2.1| ) evaluated on the trajectory ( |6.5p for j = 3/2 i.e 4 x 4 matrices, becomes 
5 = / dt(K - V) with 

= 2((4-^/3)/2 + 2/5 + 52) (6.7) 
We can diagonahse the kinetic term by using g = h — f. Then the potential term is 

V = -^ + A (2{h - 1)2 + (1 - V3)) f - 8(2 - V3)f + (10 - 4^/3)/^ (6.8) 

The plot of V{f,g) (Figs (|2|) and (^)) shows an energy landscape for a particle moving in two 
dimensions (coordinates f{t),g{t)). 



Figure 2: The potential V{f,g) for j = 3/2. 

It is clear from Fig (^) that it is possible to roll down reducible representations J, (the point 
/ = /i = 0) to the irreducible representation Jj (the point f = h = 1) without experiencing an 
energy barrier. 

In [34], we develop a geometrical picutre of this dynamical process, where two spheres merge 
into one through non-spherical deformations. 

To conclude this section, we have shown that within the ansatz ( |6.5D it is possible to roll 
down from a reducible representation to an irreducilbe representation. 

7. Effect of adding a mass term 

In this section we present briefly the effect of adding mass terms to our model ( |2.1| ): 

S = ToJ dtTii^Xf + ^[X„X,f - ^e^.kX'X^X'^ - ^X,X,) (7.1) 
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Figure 3: Contour plot ofV{f, h) for j = 3/2. 



The static equation of motion is 

[Xj, Xj] - iK €ijkXk^] - m^Xi = 0. (7.2) 

Once again this admits solutions of the form ( p. 71) 

Xi = a ©^=1 ^ (7.3) 

where now m? = 2{aK — a^), i.e. 

a = {K± \Jk^ - 2m?) /2 (7.4) 

The lower sign gives an unstable solution, so we shall only consider the upper sign henceforth. 

The model is super symmetric for m = ^ [9]. However, for the time being we will work 
with generic ra. In the table below we list the quadratic fluctuation spectrum and the value of 
the action around ([7.31) . 



Table 3: Fluctuation spectrum with mass term 



Solution 


a;2 


Multiplicity 


Action 


i©0 





7 


^(3a2 -4afi; + 3m2) 


2 

m 


3 


3(k^ + 3m^ + kVk^ - 2m?) 


8 


(h? + + K\/ K? — 2m^) 


5 


— 2m^ + K^/ K? — 2m2 


1 


2^2 





12 


4(3a2 -4aK + 3m2) 


2 

m 


9 


{k? + + K\/ K? — 2m2) 


20 


— 2m^ + K\/ — 2m2 


4 
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It is straightforward to see that we get the original energy spectrum in absence of m. On 
the other hand, when m = 2k/3 all the solutions have zero energy which is consistent with their 
being supersymmetric vacua. 

Let us rescale the field Xi / \/2 = Yi and substitute it in the energy functional of this modified 
matrix model 

^ = dt{Y^ - [Yi, Yj][Yi, Yj] + m'YiYi + h^^^.^Yi[Yj,Yk]). (7.5) 
We can write this in terms of a total square term, a total derivative term and a defect term. 

// 2\/2 i M t 4k? 

dt Y±{ -^i^Yi + -eijk[Yj. Yk\\ +TrJ dt{rr? - —)Y^ 

T Tr / dt| (^^.^ Yi + ^6,,,y,K-, n]) . (7.6) 

Lowest energy configurations are those which equate the first term to zero. This gives us the first 
order equation. Second term in the energy functional is the defect term. As it stands it is also 
a total square term but equating it to zero gives us only trivial solution. It is this term which 
vanishes in the supersymmetric limit, giving rise to the classical form of the Bogomolnyi energy 
functional. In the Bogomolnyi limit, energy functional factorises into perfect square terms and 
total derivative terms. While setting former to zero gives us first order equations of motion, 
latter gives us the energy/mass of the solution. For the time being we will treat the second term 
in the energy functional as a defect term and will not, either set it to zero or demand it as a 
constraint on the solution. The first order equations of motion are 

Yi ± (^i^Yi + '-eijk[Yj,Yk]^ = 0. (7.7) 

This equation has been studied by Bachas, Hoppc and Pioline [38]. A class of solutions studied 
by them are those, which interpolate between trivial SU(2) representation and, say, the largest 
possible irreducible representation of SU(2). The Bachas, Hoppe and Pioline ansatz which solves 
the equation of motion is 

Y^ = ^ p. (7.8) 

l + exp(V2K(t-io)/3) 

where pi is some SU(2) representation. Contribution of the total derivative term to the energy 
of this configuration comes only from the endpoints of the trajectory, i.e., from difference of the 
energy of spin representation and that of representation pi . In addition there is a contribution 
from the defect term. For the solution under consideration, this defect term contributes infinite 
energy. Recall that in the supersymmetric limit this defect term vanishes and hence there is 
no infinite contribution to the energy of this solution. Thus this is a finite energy solution in 
the supersymmetric problem. Of course, in the supersymmetric case, energy of every vacuum 
configuration is zero, and hence this solution interpolates between two degcnrate vacua. On the 
contrary, in the non-supersymmetric case, vacuum configurations have different energy which 
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depends upon their SU(2) labels. It is this non-degeneracy of these vacua which is responsible 
for the infinite energy of these solutions. For large enough A'^, even in the non-super symmetric 
scenario, we have degenerate vacua which are labelled by different set of SU(2) indices. If we 
consider a solution interpolating between these vacua, it will clearly have finite energy. 

One of the important implications of the mass term is that the moduli are lifted. This 
is evident from the first order equations. Static solutions to these equations do not allow the 
marginal deformations that we studied in section 4. Thus, it is not classically possible any more 
to roll down from a reducible representation to an irreducible one. Only way we can reach the 
irreducible representation is by tunneling. The solution discussed above falls in this class. 

8. Comments on SUGRA duals 

8.1 Polchinski-Strassler scenario 

The matrix model ( |7. 1| ) coincides with the space-independent part of the action for the ad- 
joint scalars of the four-dimensional gauge theory considered in [9]. The supersymmetric point 
m = 2k/3 and its SUGRA dual is discussed there in great detail. For other values of m, super- 
symmetry is broken; SUGRA dual of these theories also are briefly mentioned in [9]. Essentially 
the SUGRA duals of the classical solutions (|7.3| ) of the gauge theory are a collection of 5-branes 
in a space which is asymptotically AdS^ x S^. Close to themselves, the 5-brane world- volume 
is like x S'^. The radii of the 2-spheres and their location in the AdS geometry depends on 
the dimensionality of the representation jj^^^ in ( [7. 3D . 

[9] has indicated a phase diagram in the (m^ , k) plane where the < represents a 
phase in which only the irreducible representations are stable. The phase m? > 0, on the other 
hand, is qualitatively similar to the supersymmetric point m = in this phase the reducible 
representations are all stable. The matrix model (^]^) which we began with has = 0; as 
we found in our stability analysis, here the reducible representations are marginally stable, 
with classical instabilities in the cubic order. In terms of the SUGRA dual, this implies that the 
collection of 5-branes corresponding to the reducible representations should have a mild classical 
instability towards forming a large 5-brane corresponding to the irredicuble representation. 

8.2 F1-NS5 brane system 

Another place where this matrix model is relevant is in the F1-NS5 brane system in type IIA 
string theory. Easiest way to see this is to start with the D1-D5 system of type IIB string 
theory. The Chern-Simon term on the common worldvolume of the D1-D5 system contains a 
term (a, b = 0,1 are common world volume indices, i,j, k, ... = 2, 3, 4, 5 are five-brane directions, 
transverse directions ignored), 

Scs,DiD, = J rfV[. • • + $i + c%D,n^', ^^]) + ...] (8.1) 

where we have ignored terms involving higher form potentials. Since the potential C^'^^ is 
invariant under S-duality transformation, identical term exists on the common world volume of 
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F1-NS5 brane system in type IIB string theory. 



Sc.,NS = J dV[. • • + + C^]kDbn^\ ^']) + ■■■] (8.2) 

Relevant type IIA configuration of F1-NS5 branes is obtained by performing T-duality along the 
Fl string direction. Suppose Fl string is along then, T-duality along x^, leads to a term with 
clfj , which comes from C^^-j and a term with C^jl, which comes from C^^j^. Incorporating these 
effects of T-duality into (8.2), we see that the common worldvolume theory of F1-NS5 system 
in type IIA string theory contains the matrix model under consideration. The commutator 
square term comes from the DBI part of the action. To be able to do T-duality along the Fl 
direction, i.e. x^, we require all the fields (NS-NS as well as R-R) to be independent of x^. 
Therefore, effective common world volume dynamics essentially reduces down to time evolution 
only. Hence, this situation is identical to the DO brane system. 



9. Summary 

Matrix model in the RR four form background has fuzzy spheres as classical solutions demon- 
strated by Myers [4]. In this paper we studied the energy landscape of these calssical vacua. 
Contrary to our expectation we found that the reducible SU(2) representations, corresponding 
to multiple fuzzy sphere solutions are all stable. That is, quadratic fluctuations around these 
solutions do not have tachyonic instability although these solutions all have higher energy the 
the irreducible representation. We resolve this puzzle by showing that the equations of motion 
allows many more solutions. These new solutions are energetically degenerate with the known 
solutions, uu. SU(2) representations. In other words, original solutions have a large moduli 
space. Quadratic fluctuations around the new deformed solutions do have tachyonic instability 
indicating the roll down path towards more stable configurations. The irreducible representation 
of SU(2) has no nontrivial flat direction as well as the quadratic fluctuation does not have any 
tachyonic modes. This feature is expected as the irreducible representation is the lowest energy 
solution. We then compared our results with other approaches, particularly the spherical D2 
brane of Bachas, Douglas and Schweigert [6], WZW matrix model of Alekseev, Recknagel and 
Schomerus [12], and Myers dual D2 brane in the presence of RR four form flux. 

We then discussed the effect of marginal deformation and studied the dynamical trajectory 
which showed that it is classically possible to roll down from a configuration corresponding to a 
reducible representation to an irreducible representation. We also studied the effect of switching 
on the mass term in the matrix model. For a specific value of this mass perturbation, the matrix 
model becomes super symmetric. At this point all the classical vacua become degenerate with 
zero energy. Precisely at this point the second order differential equation of motion can be 
factored into a first order equation. The main new observation in this case is that the moduli 
are lifted due to mass perturbation. Thus, it is not classically possible any more to roll down 
from a reducible representation to an irreducible one. However, time dependent solutions to the 
equation of motion do exist. These tunneling solutions non-perturbatively interpolate between 
different vacua. 
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There are several leads which are worth pursuing. Firstly, full time-development of the roll 
down solution in the massless case is important particularly from the viewpoint of cosmology. 
Along the same lines it will be interesting to develop the time-dependent solution in supergravity 
(cf. Sec. 8). Secondly, there is an interesting issue related to topology change due to non- 
commutativity. In the massless case, there is a topology change classically [34]. On the other 
hand, in the massive case, there is no classical path interpolating between different classical 
vacua and therefore, there is no topology change classically, but it can occur via quantum 
tunnelling [39]. 

Acknowledgement: We would like to acknowledge discussions with Sumit Das, Shiraz Min- 
walla and Sandip Trivedi. 



21 



A. Eigenvalues 

In this appendix, we calculate the eigenvalues 
K^uP' of the quadratic fluctuation matrix M 
in ( |3.1| ) around the general, static classical so- 
lution ([4.1|). We present the calculations for 
n X n matrices for n = 2, 3, 4, 5. 



N=2 


Solution 




Multiplicity 


0©0 





5 




4(±c + c2) 


2 (each) 


1/2 





3 




2 


6 



Here c is the separation between the zero branes, 
i.e the solution is Xi = Cia^. 



N=3 


Solution 


2 


Multiplicity 


oeo® 





12 


4(±c + c2) 


2 (each) 


\c-d^±2\c^d 


2 (each) 


\c + d^ ±2\c + d\ 


2 (each) 







10 


2 


6 


|(l + 4|cl2±4V3|cl) 


2 (each) 


i(3 + l2|cT±4V3|c1) 


2 (each) 


1 





8 


2 


6 


6 


10 



The solution in the 3 0-branes case, is defined 
as follows: 

Xi = CiXs + ^Ag, 
where the As^g's are the Cartan elements. 
In the spin 1/2 + 0-brane case, the solution is 

Xi = Xi + CiXs, 
where the Aj are non-zero only in the upper 
left 2x2 block, and form an su(2) sub-algebra. 



N=4 



Solution 


, ,2 
OJ 


iViultiplicity 




U 


1 o 




o 
Z 





2 © U ® 


j[l ± 4Voc + 4c ) 


z (eacnj 




^(d ± 4v oc + 12c j 


z (eacnj 




4 /'II ^11 2 I /oil ^in 


z (eacnj 




62 


z (eacnj 




U 


1 7 
1 i 




o 
Z 


Iz 


2 2 


4(±c + c ) 


z (eacnj 




z(i ± 4c + zc j 


z (eacnj 




Z[i ± zc + Zc ) 


z (eacnj 




2(1 + 2c j 


4 




U 


1 

1 ( 




2 


D 









1 ® u 


■^[6 ± zyoc + 4c j 


z (eacnj 




2 / o _|_ /I , _i_ /I ^2 \ 
^ t^O =t 4 V DC -f- 4C j 


z (eacn j 




1(3 + 4c^) 


4 
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3 
2 


2 


6 




6 


10 




12 


14 



\\af = \^'^ + 2\d\'^-2^/2c■d = +d'^ -2V2cd. 

62 = (1/12)(9 + 4(c + 2^/2d)2 ± 4^/3(c + 
2^/2(i) ±8V3(c + 2^/2d)). 
The solution in the spin 1/2 + 2 0-branes case, 
is defined as follows: Xi = Xi + cXs + dXi^, 
where Xi are the spin- 1/2 rotation generators 
in the upper left 2x2 block (and zero in the 
rest) and As and A15 are the Cartan generators 
which are equal to identity in the upper left 
2x2 block. This situation represents an 5^ 
brane at the origin and two DO-branes located 
at positions labelled by c and d. 
In the spin 1/2 + 1/2 case, the two spin half 
2-branes are separated along their center of 
mass by the vector c, i.e the solution is 

Xi = ai (g) I + Cil (g) (T3 
Similarly, in the spin 1+0— brane, the solution 
is 

Xi = Ji + Cj Ai5 
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where Ji are the spin-1 rotation generators in the upper left 3x3 block, and A15 is the Cartan 
element (l/-v/6)L>m(7(l, 1, 1, — 3). 



N=5 


Solution 




Multiplicity 


i ® 0® 





42 


ffi 


3/4 


24 




2 


6 







28 




z 


D 







lU 




(z/3j(3 + 4c j 


4 




(_Z/0j(_0 ± ZV DC + 4c j 


z (^eacnj 


1 U (13 U 


(_Z/0j(^0 ± 4V DC + 4c j 


z (^eacn j 




(1/6)(12+ ) 


4 




(l/6)(3|7l ± 2V6||7lr) 


2 (each) 




(1/6)(12±2V6|/3| + \I3Y) 


2 (each) 




(1/6)(12±4V6|/3| + \I3Y) 


2 (each) 







32 




3/4 


16 




2 


24 







26 




3/4 


8 


i®i 


2 


12 




15/4 


16 




6 


10 



||^||2 = ||^|2 + (5/3)||rf||2_2y'|c.d 

\\Pf = \\4'' + lf>\\df + 2Vl^c-d 
and c = ||c||. 

Note that the solution in the spin 1 + 2 0-branes case, is defined as follows: Xi = Aj + cA2i + (iA24, 
where Aj are the spin one rotation generators in the upper left 3X3 block (and zero in the rest) 
and A21 and A24 are the Cartan generators which are equal to identity in the upper left 3X3 
block. This situation represents an S"^ brane at the origin and two DO-branes located at positions 
labelled by c and d. 
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